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Abstract 

The experimental observation of effects due to Berry's phase in quantum systems 
is certainly one of the most impressive demonstrations of the correctness of the 
superposition principle in quantum mechanics. Since Berry's original paper in 1984, 
the spin 1/2 coupled with rotating external magnetic field has been one of the 
most studied models where those phases appear. We also consider a special case 
of this soluble model. A detailed analysis of the coupled differential equations 
and comparison with exact results teach us why the usual procedure (of neglecting 
nondiagonal terms) is mathematically sound. 
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The study of quantum systems in the adiabatic regime appeared for the first time 
in the paper by Born and Fock in 19280. Since their original paper it is known that 
the time evolution of an instantaneous eigenstate of the time- dependent Hamiltonian 
acquires a phase besides the dynamical phase. But it was only after Berry's seminal 
paper in 1984 || that this extra phase was raised to the condition of a geometric phase. 
Different quantum systems in this regime have attracted a lot of attention due to the fact 
that the existence of those geometric phases introduces a shift in the frequency initially 
predicted and afterwards experimentally measured. There is a list of experimental articles 
on observation of Berry's phase in Ref . || . More recently, geometric phases were observed 
in neutrons 0] and in an atomic state interacting with a laser field ||. 

Even though the research on quantum systems in the adiabatic regime is still a very 
active area in physics nowadays, these geometric phases have become an issue included 
in undergraduate textbooks in quantum mechanics^]. Certainly this comes from the fact 
that the measured frequency shift due to the presence of geometric phases is one of the 
proofs of the linearity of the quantum phenomena. 

The traditional characterization of a quantum system in the adiabatic regime is that 

its interaction with the environment is described through a set of classical time-dependent 

parameters R(t) = (Ri(t), ■ ■ ■ , R n (t)) periodic in time but whose period is much larger 

than the characteristic time scale involved in the quantum phenomena. The Hamiltonian 

that drives the quantum system depends on the set R(£), that is H[R(i)]. Let bj be 

the angular frequency of R(£) and oji the angular frequencies equal to the difference 

between two distinct quantum eigenvalues of H[R(t)] divided by h. The adiabatic regime 

is attained when ^ — > 0, for all a;, 0, || [7] . Under these conditions, we have the so-called 

the Adiabatic Theorem^: 

For a slowly varing Hamiltonian, the instantaneous eigenstates of the Hamiltonian 
evolve continously into the corresponding eigenstate at a later time. 

Let us consider a periodic Hamiltonian that returns to its original operator, that is, 
H[R(T)] = H[R(0)] when the external classical set R(£) completes a period. If the initial 
quantum state is an eigenstate of H[R(0)], the adiabatic theorem affirms that the evolved 
state is an eigenstate of the instantaneous Hamiltonian with the same initial quantum 
number but it is does not forbid it to acquire a phase. This phase can be decomposed 
into a dynamical phase plus another one that M.V. Berry proved to be geometric f2j. This 
geometrical phase acquired by instantaneous eigenstates of the Hamiltonian chosen as 
initial state of the quantum system is more easily recognized if the time evolution of the 
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initial state is described in the basis of the instantaneous eigenstates of the Hamiltonian 
H[R(£)]. To make our notation simpler, we denote: H[R(t)] = H(£). Let 1^(0)) be the 
initial state of the quantum system; its time evolution is given by \ifi(t)). We decompose 
\ip(t)) on the basis of the instantaneous eigenstates of the Hamiltonian: 



N ■ rt 

m)) = Y, c n{t)e- iS ° En{t)dt ' \M, (1) 



n=l 



where \4> n ',t) are the instantaneous eigenstates of H(t) with eigenvalue E n (t), 

H(t)\ ( f )n ,t) = E n (t)\ ( p n ,t). (2) 

The index n represents a set of quantum numbers that uniquely determine this quan- 
tum state. We are assuming that the time-dependent Hamiltonian has iV nondegenerate 
instantaneous eigenstates. The phase f* E n (t')dt' is called the dynamical phase and 
reduces to ~ lEnt when the eigenvalues of the Hamiltonian are time independent. Requiring 
expansion ([[]) to be a solution to the Schrodinger equation 

m\m)=in^-, (3) 

we get a set of coupled differential equations for the coefficients c n (£)[@, ^, [jj: 



^ = -(M^tMt) + £ ^'^e^ eU^- E ^' c m (t). (4) 

From Eq. QT|) up to Eq.(f|) the expressions are exact. In Refs. [El, [7|, the adiabatic 
approximation is implemented by recognizing that 

I A +1 9 U +\ (0n;*lnw^l0m;*) f , /r\ 

Wn;t\ — \(j) m ;t} = at , for n^m. (5) 

at £j n — iL m 

In this approximation those terms are neglected compared to the elements (<f) n ; t\^\(j) n \ t) 
in Eq.(||). After neglecting those terms, the differential equations (f|) decouple and each 
one is easily solved || 

However, since the non-neglected terms in Eq.(^) are also a measure of how much the 
instantaneous eigenstates of the Hamiltonian vary with time, they must also be small once 
this variation in time is due to the interaction with the environment and by hypothesis 
it is very slow. At first sight this assertion would seem to contradict the traditional 
derivation of the results in the adiabatic limit p], |[ [7|]. The aim of this note is to show 
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that this does not happen. We explicitly show this in the context of a soluble model. 
In the differential equations for the coefficients of the expansion of the vector state [see 
Eq. ([!])] for this model we introduce tracers to follow the contribution of crossed terms 
to the exact solutions and finally we explicitly show why the adiabatic approximation is 
correctly obtained by only neglecting the terms (|5]). 

We consider a particular case of the model presented by M.V. Berry in Ref. @, that 
is, a spin | coupled to an external magnetic field B{t) with constant norm B (B = |J3(£)|) 
that precesses with constant angular frequency u> around the z-axis: 

B(t) = (Bsin(9)cos(u t), Bsin(9)sin(uj t), Bcos(9)) . (6) 

All directions in space are equivalent, and we assume that < 9 < ~. 
The Hamiltonian of the model is 

H(t) = /iB(t) ■ s. (7) 

The constant /i is written in terms of Lande factor g and Bohr magneton /iefl, that is, 
li = g^B and s is the spin | operator. Hamiltonian (0) is the interaction energy between 
the intrinsic magnetic moment of the electron and the external magnetic field0. Its 
instantaneous eigenstates are: 



n;t) = -A\ t) + De lu)ot \ [) Ex 
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\<p 2 -t) = D\])+Ae^ t \{) =► E 2 = ^, (9) 

where | |) (| J,)) is the eigenstate of the component s z of the spin operator with eigenvalue 
| (— |). The state \(f>i;t) (\(f)2', t)) is the eigenvector of spin down (up) along the direction 
n = ^jp-. We defined 

9 9 
A = sin{-) and D = cos(-). (10) 

The eigenvalues E\ and E 2 are time independent. For this model, Eq.(0) becomes 

\m) = jlc n {t)e-^ t o E ^ )dt '\ ( t> n} t). (11) 

n=l 

The coefficients c\(t) and Oz{t) satisfy the differential equations 

= -iu D 2 Cl {t) - tu ADe- 2i " lt c 2 {t) (12) 

LLC 
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= -iu Q ADe 2luJlt Cl {t) - iu A 2 c 2 (t). (13) 

In this model we have only one characteristic frequency of the quantum phenomena. It 
is proportional to the difference of the eigenvalues of Hamiltonian (|7|), 

2wi = - = (14) 

The factor 2 in definition (|14]) is introduced for convenience. In Eq.(|I~2"D the adiabatic 
approximation (uq — > 0) implies that the C2(t) on its rhs has to be neglected. However, 
the coefficient that multiplies c±(t) on the rhs of the same equation is also proportional to 
lvq although it is not neglected. A similar discussion is valid for Eq. ([13|) by interchanging 
C\(t) and c 2 (t). We should remember that in the adiabatic approximation we have no 
restriction on the #-angle, except that it has to be different from 9 = 0° and 9 = n. 
At first glance it seems that we are not consistenly keeping the terms in our coupled 
differential equations ( |T2"D and flT3]). We will see that is not the case and the difference 
relies on the exponential that multiplies 02(f) [ci(t)] on the rhs of Eq.fll2|) [Eq.(|K 



To follow the contribution due to each term on the rhs of Eqs. ([12]) and fllcf) , we 
introduce tracers for each term on the rhs in the equations, that is, 

dCl ^ ■- -ianuJoD 2 ^) - iau ADe- 2lu)lt c 2 {t) (15) 



dt 

-iaLO Q ADe llUlt ci{t) - ia 22 uJoA 2 c 2 {t) , (16) 



dt 

where an, 022 and a are the tracers. These coupled equations are exactly solved. At the 
end we take the tracers equal to 1 or depending on the approximation under considera- 
tion. The differential equations in the adiabatic approximation are recovered by chosing 
a = 0. To make the rhs of the equations time independent, we define new variables, which 
is equivalent to going to the rotating frame. The new variables are 

Xi(t) = ci^y* 1 * and X 2 (t) = c 2 (t)e- iu;it , (17) 
and they satisfy the equations 

= i( UJl -a u u D 2 )X 1 (t)-iaoj ADX 2 (t) (18) 

dt 

= -iau) Q ADX 1 {t)-i{uj 1 + a22U)QA 2 )X 2 (t). (19) 



dt 
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Taking the general initial condition 



|^O)> = Cl (O)|0 i; O) + C2 (O)|0 2 ;O>, (20) 

where |</>jj 0), i — 1, 2, are the eigenstates of Hamiltonian (^) at t — 0. The exact solutions 
of Eqs. (p~8|) and (0) with the initial state (|(J) are 



d(t) 



-i[^l + ^l(aiiD 2 +a22A 2 )]t 



{cx(0)cos(Tt)+ 



and 



+ 



z sm(rt) 



(u x + ^-(a 22 A 2 - a nJ D 2 )) Cl (0) - au ADc 2 (0) 



(21) 



c 2 (t) 



Ji~^(a 11 D 2 +a 2 2A 2 



{c 2 (0)cos(Tt) + 



and 



z sm(rt) 



(u x + y (a 22 A 2 - a n D 2 ))c 2 (0) + aio AD Cl (0) 



(22) 



r = lji 



l + ^(a 22 A 2 



_, 1 

_> 



(23) 



From Eqs. (|2l|) and ([22]) we see that the nondiagonal terms in Eqs. (|15|) and (|16|) (a ^ 0) 
do not contribute to the overall phase in the expressions of c\(t) and c 2 (t). On the other 
hand, in the expressions of these coefficients there appears a new frequency T (Rabi's 
frequency |10"|, [□]]). For u>\ ^ and u>\ lu , the contribution of the diagonal terms 
(an 7^ and a 22 ^ 0) to the expression of frequency T is O(^) while the contribution 
of the nondiagonal terms (a ^ 0) is 0({^) 2 ). We get the geometric phase0 when the 
external parameter (in our case B(t)) completes a period t — ^0 and for this time 
scale we can neglect contributions of 0((^) 2 ) to the frequency T. Those are exactly the 
contributions coming from nondiagonal terms in Eqs. (p~5| ) and fli~6|). Due to this fact, we 
get the adiabatic limit by substituting a = in the expression of T. 

The nondiagonal terms in Eqs.flTo]) and flTS| ) also contribute to the sine-term in Eqs. 
(f2"T|) and (|2"2"|). We now discuss the sine-term on the rhs of Eq.fl2~T|). A similar discussion 
is valid for Eq . (^2|) . The sine-term in Eq.(21) is rewritten as 



isiniTt) 



l-a\^fA 2 D 2 yci(0)-a^ADc 2 (0) 



(24) 
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In the adiabatic regime we have ^ < 1 and then from Eq.(p3|) it follows that ^ < 1. 
Keeping terms up to zeroth order in ^ 1 in the sine-term (123), we obtain 



isiniTt) 



l-a 2 (^) 2 A 2 D 2 y Cl (0) - a^ADc 2 (0) 



Substituting the approximation 
limit 



ici(0)sin(rt). (25) 
in the expression of c\(t), we get in the adiabatic 



c"- a -(t) = exp 



-i ( ui + — [anD + (I22A 



with 



l + ^{a 22 A 2 -a n D 2 



that is, the Rabi's frequency in the adiabatic approximation. 
By the same procedure we obtain, in the adiabatic regime, 



ci(0), 



(26) 



(27) 



c 2 a -(t) = exp 



l [ UJl -^-(a u D 2 + a 22 A 2 )-r a - a -)t 



c 2 (0). 



:28i 



Taking an = a 22 = 1, in Eqs.(p6"D and (^) we recover the known results for the adiabatic 
approximation |2|, ||, which is equivalent to making a = in Eqs.(|15D and (|IIJ) from the 
beginning. 

We showed analitically, for a two-level model, why we get the correct adiabatic limit 
expressions by just neglecting the nondiagonal terms in the coupled equations ( |P2"|) and 
(|T3|). These first-order differential equations drive the time evolution of the expansion 
coefficients flip for the quantum state on the basis of the instantaneous eigenstate of the 
time-dependent Hamiltonian. 

It certainly could be questioned why terms of different orders in (^) are dropped in 
the expansions of the coefficients Cj(£), i = 1,2, and of the frequency T. The big difference 
between these two expansions is due to the fact that we are interested in time scale of 
order T — — , where the linear correction in (— ) to T gives a finite contribution. This 
does not happen to the linear correction (^ a ) to the amplitude of the terms 



sin(Tt) 



Even though we studied a soluble spin | model, many features of this model appear 
in any quantum system that is coupled to a tridimensional external parameter that varies 
in time very slowly. For these models we can have a matrix representation of Eq.(^) on a 
suitable basis such that the coefficients can be written as a column where the equations 
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that drive their time evolution are a generalization of Eqs.(|i~2|) and (jlB]). Neglecting the 
nondiagonal terms in the adiabatic limit is justified if they are multiplied by a phase whose 
angular velocity is much larger than the angular frequency of the external parameter. 
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